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Abstract Improvement of the simple tomography method 
for estimating moisture distribution in squared timbers was 
conducted by using impedance measurement in the RF 
range at 250 kHz. The local impedance in the horizontal and 
vertical directions was measured during natural drying. An 
inverse procedure using a trial-and-error method was used 
to reconstitute the cross-sectional moisture distribution 
from the surface impedance measurement. Two essential 
corrections for the spatial heterogeneity of the electric field 
between electrodes and for the calculation error due to the 
algorithm of the trial-and-error method were introduced to 
improve the accuracy of the reconstitution of the moisture 
distribution. The reconstituted moisture distribution 
matched well that determined by the oven-dry method. 
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Introduction 

Many types of electronic moisture meters are used in the 
wood industry. However, these meters are not sufficient to 
detect the moisture gradient in structural timbers. A com¬ 
puter tomography method using X-rays has been studied 
for standing trees * 1 and for square timbers. 2 Another method 
using neutron beam imaging has been studied for wood 
disks. 3 These techniques can reproduce the cross-sectional 
moisture distribution, but they are too expensive to intro¬ 
duce into the wood industry. The other tomography method 
using impedance measurement has been studied in the high- 
frequency range 4 8 and using direct current. 9 From the point 
of view of improvement of the conventional electronic 
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meter, a pin-type sensing probe with multiple hollow elec¬ 
trodes indenting in a specimen has been designed for mea- 
suring the moisture gradient at different depths. Our final 
target is application of the impedance tomography tech¬ 
nique in the radio frequency range to the inline grading of 
structural timbers with a low-cost moisture meter. 

A recent study 4 showed potential for estimating the 
cross-sectional moisture distribution in structural timber by 
using the impedance tomography method. However, the 
reconstitution of the moisture distribution was not ade¬ 
quate for practical commercial use. This work aims to 
improve the reconstitution of moisture distribution by con¬ 
sidering two kinds of signal correction. One corrects for the 
effect of the heterogeneity of the electric field between 
the electrodes and the other corrects for the flattening of 
the peak of the distribution pattern due to the inverse cal¬ 
culation algorithm used. 


Theory 

Reconstitution of moisture content distribution 

The algorithm used for the inverse calculation is commonly 
used for inverse problems in acoustic tomography and has 
the benefit of simple calculation but with limited resolu- 
tion. In this article, the moisture distribution was deter¬ 
mined so that the moisture content of individual cells in Fig. 
1 satisfy Eq. 1 using the simple trial-and-error method. 
Figure 1 shows the cross data table of 6 x 6 elements used 
to obtain the inverse solution. The value in each cell u if is 
the local moisture content. U Hj and U yi are accumulated 
moisture contents of the rows and the columns, 
respectively: 

and (A,=i>„ (1) 

7=1 i =1 

The procedure for obtaining the solution is as follows: 

In the first step, an average value of each accumulated 
moisture content, U Ui , is distributed to each row cell: 
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Fig. 1. Cross-sectional data table of 6 x 6 elements of moisture content 
distribution. U Hi and U Vj are the accumulated moisture contents in the 
row and column directions, respectively 


Detection side electrodes 



Fig. 2. Arrangement of the electrodes. GE, guard electrode 


u’i=U w /6 (2) 

In the second step, the average of the difference of the 
accumulated moisture content in the column direction 
between U Vj and the first estimated value is distributed to 
each row cell: 

u 

This procedure is repeated in the row direction, and these 
two procedures are iterated until each value in the cross 
table converges to a stationary value. These series of inverse 
calculations were executed by an Excel (Microsoft) macro 
program. 


"=^ + fc/ v/ -X^ 

V i =1 J 



( 3 ) 


Materials and methods 

Boxed-heart squared green hemlock (Tsuga heterophylla) 
timber of 105 mm x 105 mm x 4 m was used. It was cut into 
twelve short specimens of 320 mm and both ends were 
sealed with silicone sealant. The specimens were dried 
under room conditions in the temperature range 13°-20°C 
and the average air dried density of the specimen was 
490 kg/m 3 . An impedance meter (Hioki-3531 Z HiTESTER) 
was used. The absolute value of the wood impedance Z was 
measured at 250 kHz with an excitation voltage of 1 V and 
the electrode arrangement shown in Fig. 2. The local imped¬ 
ance was measured by switching the detection electrode in 
turn from positions 1 to 6. A pair of thin polyethylene sheets 
was sandwiched between the electrodes and the timber. 
After the electrical measurement, the specimen was cut into 
36 (6 x 6) pieces to determine the moisture content distribu¬ 
tion by an oven-dry method. 
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Fig. 3. Relationship between impedance ln(Z) and accumulated mois¬ 
ture content {MC) U of pooled data. Open circles, row direction; solid 
circles, column direction; R 2 , coefficient of determination; solid line, 
column direction; dashed line, row direction 


Figure 3 shows the relationship between the impedance 
and the accumulated moisture content in the directions of 
the 6 rows and 6 columns for the pooled data through all 
drying stages, where the accumulated moisture content is 
defined as the sum of moisture contents in each row or 
column direction. There was no significant difference 
between the regression equations for the row and column 
directions. Since these equations have acceptable correla¬ 
tion coefficients, the data in both row and column directions 
are pooled in the following analysis. 


Corrections improving reconstitution of 
moisture distribution 


Results and discussion 

Fundamental equation converting impedance to 
moisture content 

Past work 4 showed that the impedance measurement at 
250 kHz is appropriate for the moisture measurement 
because the impedance change due to the moisture content 
change is stable and large at this frequency. 


It is well known and has been experimentally demonstrated 
that the temperature rise of timber piled between elec¬ 
trodes used for high-frequency drying of wood depends on 
the spatial position. 12 Hence the strength of the electric field 
between the electrodes shows a spatial heterogeneity. 

Accordingly, prior to the reconstitution of the moisture 
distribution from the impedance data, the following 
two corrections were carried out: first, the correction for 
heterogeneity of the electric field in the cross section of the 
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Fig. 4. Effect of spatial location in the relationship between impedance 
ln(Z) and U. Solid circles , core group; open circles , outer group 

material, and second, the correction for the inverse 
calculation. 

Correction for heterogeneity of electric field 

The first correction is to compensate for the field heteroge¬ 
neity and it would depend on the spatial configuration of 
the electrodes, specimen size, and discontinuity of material 
properties. For the particular geometrical conditions of this 
experiment, the following two spatial groups were consid¬ 
ered. One is made up of the outer rows and outer columns 
(£/ vl , U y6 , U m , and U m in Fig. 1), and the other is made up 
of the core rows and columns (U V2 -U V5 and U m -U H5 in Fig. 
1). The impedance versus moisture relationships were 
obtained for these two distinct groups. 

Figure 4 shows the relationship between the impedance 
Z and the accumulated moisture content U. The regression 
lines for the core group and outer group are as follows: 

Core group: ln(Z core ) =-0.0100C/ +14.112 R 2 =0.635 (4) 

Outer group: ln(Z out ) =-0.0151C/ +14.394 R 2 = 0.760 (5) 

where R 2 is the coefficient of determination. The correlation 
coefficients were improved over those in Fig. 3 by separat¬ 
ing the data into the two groups. 

The estimated values of the accumulated moisture 
content Z ESX in the core and outer parts are given by the 
following equations inversely from Eqs. 4 and 5, 
respectively: 

^ESTcore = 1411.2-1001n(Z core ) (6) 

^E STout = 953.2-66.231n(Z out ) (7) 

Figure 5 shows the relationship between the estimated accu¬ 
mulated moisture content U EST for both row and column 
directions and the oven-dry accumulated moisture content 
U oven? and the regression equation is given by: 

t/ OVEN = 0.7522 £/ est + 49.757 R 2 = 0.752 (8) 

Although the equation fits the data reasonably well ( R 2 = 
0.752), the predicted moisture content is under-estimated 
when below 200% and over-estimated when above 200%. 



Estimated AMC (%) 

Fig. 5. Relationship between oven-dry accumulated moisture content 
U oven and estimated accumulated moisture content U EST . AMC, accu¬ 
mulated moisture content 


In order to correct this trend, the combination of the Eqs. 
6-8 was altered as follows: 

t/ESTcore =1111.3-75.221n(Z core ) R 2 =0.7523 (9) 

f/ESTout =766.8-49.811n(Z out ) R 2 =0.7523 (10) 

This treatment, while unorthodox for linear regression, 
addressed the nonlinearity inherent in the electric field dis¬ 
tribution. The resulting relationship between the modified 
moisture content and the oven-dry moisture content was 
given by: 

XJ OVEN = XJ £ ST — 0.0015 R 2 = 0.7523 (11) 

The gradient of the regression line became 1, and there 
was no difference in the coefficient of determination from 
Eq. 8. 

Correction for algorithm in inverse calculation 

The second correction in the inverse calculation was also a 
development beyond the previous work, 4 where the itera¬ 
tive method used was apt to flatten the peak of the resolved 
distribution. The effect was quantitatively evaluated using 
the original data set obtained from the cross-sectional dis¬ 
tribution of the oven-dry moisture content. The original 
oven-dry moisture content data were compared with the 
reconstituted values to obtain the correction factors from 
eight data sets. The results were averaged and are shown in 
Fig. 6, which gives the magnitude of the correction factor 
applied to sharpen the moisture distribution. The correction 
factor in Fig. 6 shows that the uncorrected iteration method 
has a tendency to underestimate the core and four corner 
sections and to overestimate the remainder. 

Application of two corrections to impedance data 

The two corrections described above were adopted in the 
following order: first, the corrections for spatial position 
given by Eqs. 9 and 10 were applied and then the correction 
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for the inverse calculation procedure using the correction 
factor in Fig. 6 was adopted. Figure 7 shows the moisture 
content distributions for the oven-dry method and for the 
tomography method without corrections and after correc¬ 
tions. The prediction of moisture content distribution using 
the corrected inverse calculation showed a much closer 
alignment with the oven-dry results than those without the 
corrections. 

Figure 8 shows the moisture content distribution for the 
oven-dry method and the reconstitution method. The rela- 
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tionship of the local moisture contents between the oven- 
dry method and the reconstitution method was given by the 
following equation: 

Woven = 1-005w ESX - 0.582 R 2 = 0.713 (12) 

The deviation of Eq. 12 from the oven-dry results became 
large above the fiber saturation point, but since the main 
object of the application is the moisture inspection of struc¬ 
tural timbers below the fiber saturation point, the obtained 
results are sufficient for practical evaluation of the moisture 
gradient of commercial timbers. The other important point 
is that the reconstituted cross-sectional moisture distribu¬ 
tion gives us useful information to evaluate the drying 
process by determining whether the timber moisture content 
is within the target range to assure product quality. If a large 
moisture gradient remains in the timber, it may induce deg¬ 
radation such as checks and crooks during service. Figure 9 
shows the relationship between the oven-dry average mois¬ 
ture content in the whole cross section and the predicted 
moisture content. The regression equation is as follows: 

U whod = 0.8678t/ whest + 4.780 R 2 = 0.780 (13) 

An alternative regression line with zero intercept is as 
follows: 

U whod = 1.007f/ whest R 2 = 0.759. (14) 


Fig. 6. Correction factor per unit moisture content applied prior to 
inverse calculation 


The average moisture content is also well estimated by this 
tomography method. 
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Fig. 7. Reconstituted moisture content (%) distributions. Left, oven-dry moisture content; middle, no correction; right, after both corrections 
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Fig. 8. Relationship between w OVE n and w EST for individual positions 
during all drying stages 



10 20 30 40 50 

Estimated whole MC (%) 

Fig. 9. Correlation of moisture content of whole specimen between the 
oven-dry method I/ whod and estimation £/ whest 

Further improvement to the described method would 
require another calculation algorithm that takes account of 
the curvilinear nature of the electric force lines. Introduc¬ 
tion of information on diagonal elements such as by using 
X-ray computer tomography methods could be used, but 
these are complex and expensive. The advantages of the 
impedance measurement are not only low equipment and 
processing costs, but also safety and ease of handling in the 
field and in the saw mill. 


Conclusions 

The proposed method aimed at improvement of the recon¬ 
stitution of the moisture distribution in timber using imped¬ 


ance measurements. Corrections for the effects of the 
heterogeneity of the electric field and the inverse calcula¬ 
tion algorithm were introduced. The results showed better 
reconstitution of the local moisture contents of commercial 
timbers than the previous approach did. This method is a 
key technique leading to realization of a moisture meter 
estimating the cross-sectional moisture distribution in com¬ 
mercial timbers. 
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